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Overview 


The research carried out during this grant achieved the scientific goals as set forth in the initial 
proposal. By way of introduction, this research project was initially awarded by NASA to the Principal 
Investigator (PI) when he was employed by The Aerospace Corporation; the project also has a co-Principal 
Investigator, Dr. Ramon Lopez, whose funding was separate. The initial proposal was a three-year effort, 
the first year of which was completed while the PI was an employee of The Aerospace Corporation. During 
the last year s(coincident with the start date of the period covered in this report), the PI changed institutional 
affiliation from The Aerospace Corporation to Boston University (BU). Following the move, the PI 
continued the work on the overall research effort described herein by establishing a “new”, restructured 
two-year grant at BU whose tasks included those of the remaining two years of the initially proposed study, 
and which included a no-cost extension. This final report documents primarily those elements not 
described in previous progress reports. 

PROJECT HlSTORY/EVOLimON 

To put the results of this effort into perspective, it is first important to provide the project context and 
history owing to the institutional changes and task restructuring. The main thrusts of the overall proposed 
effort are to understand the three-dimensional structure of the substorm current wedge, its magnetospheric 
sources, and the magnetic field geometry throughout the magnetotail during the substorm cycle. Our 
approach has been one of progress through iterative use of data analysis, modeling, and theory. In the first 
year of the effort, we focused on three main tasks (principal collaborators are shown parenthetically) : 
establishing the magnetospheric equatorial projection of the ionospheric electrojets using empirical, static 
magnetic field models (Dr. Ramon Lopez) ; building an improved, pre-onset magnetic field model of the 
ring current and near-magnetotail region by appealing to empirical, equilibrium plasma pressure profiles 
(Dr. Anthony T. Lui, Dr. David P. Stem); and including explicitly the effects of field-aligned currents in 
magnetic field models in order to assess their importance when mapping ionospheric currents to their 
magnetospheric source regions (Dr. Eric Donovan, Dr. Ramon Lopez). Several presentations of these 
multi-faceted efforts were reported at scientific meetings and were completed and published. 



Task Transition 

During the final period of the effort, all principal collaborators (as well as the PI) of the outstanding 
tasks experienced institutional changes (Donovan, Lopez, and Spence) and the associated shifts in priorities 
and duties. This wholescale reorganization led to logistical problems of focusing jointly on the specific 
remaining tasks begun in the first year that had originally been anticipated for year two of the effort We 
reassessed our positions and decided to focus more tightly on the individual pieces of our study during this 
final part of the funding cycle with the intent that the conclusion of the outstanding tasks begun in the first 
year of the effort be continued at a low level and finished during the final year of the grant 

Final Summary 

This section outlines those tasks undertaken in the final year that contribute integrally to the overarching 
project goals. First, during the final year, it is important to note that the project benefited greatly with the 
addition of a Boston University graduate student, Ms. Karen Hirsch. Jointly, we made substantial progress 
on the development of and improvements to magnetotail magnetic field and plasma models. The ultimate 
aim of this specific task was to assess critically the utility of such models for mapping low-altitude 
phenomena into the magnetotail (and vice-versa). The bulk of this effort center ed around the finite-width- 
magnetotail convection model developed by and described by Spence and Kivelson (7. Geophys. Res., 98, 
15,487, 1993). This analytic, theoretical model specifies the bulk plasma characteristics of the magnetotail 
plasma sheet (number density, temperature, pressure) across the full width of the tail from the inner edge of 
the plasma sheet to lunar distances. Model outputs are specified by boundary conditions of the source 
particle populations as well as the magnetic and electric field configuration. 

During the reporting period, we modified this code such that it can be interfaced with the auroral particle 
precipitation model developed by Dr. Terry Onsager. Together, our models provide a simple analytic 
specification of the equatorial distribution of fields and plasma along with their low-altitude consequences. 
Specifically, we have built a simple, yet powerful tool which allows us to indirectly “map” auroral 
precipitation signatures (VDIS, inverted-V’s, etc.) measured by polar orbiting spacecraft in the ionosphere, 
to the magnetospheric equatorial plane. The combined models allow us to associate latitudinal gradients 
measured in the ion energy fluxes at low-altitudes with the large-scale pressure gradients in the equatorial 
plane. Given this global, quasi-static association, we can then make fairly strong statements regarding the 
location of discrete features in the context of the global picture. We reported on our initial study at national 
and international meetings and published the results of our predictions of the low-altitude signatures of the 
plasma sheet 

In addition, the PI was invited to contribute a publication to the so-called “Great Debate in Space 
Physics” series that is a feature of EOS. The topic was on the nature of magnetospheric substorms. 

Specific questions of the when and where a substorm occurs and the connection between the auroral and 
magnetospheric components were discussed in that paper. This paper therefore was derived exclusively 
from the research supported by this grant 
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Introducti on 

Although the magnetosphere has been studied and sur- 
veyed extensively for over 3 decades since the advent of 
th6re are ^outstanding tasks pertai^ £ 
U “ e magnetosphere. One of these is aquietri me 
magnetosphenc model with sjxxification of the magneric 
field configuration and the associated equilibrium pEma 
pressure distribution. In the quiet time Earth's magneto- 
sphere, the inertial force pdv/df is usually insignificant 
d therefore force balance should exist between magneto-’ 
spheric plasma pressure and electromagnetic fon^A^ 
explicit configuration resembling the mae-netn^nhor* a 
satisfying force balance has yet to £ IS 

Equilibrium of the magnetosphere has been studied 
theoretically in the two limiting situations of the inner 
magnetosphere, where the Earth’s dipole field doming 

Shi iK ! where th * field arises mainly 

the theta-shape tail current system. The inclusion^ the 

transition region between these two simplified situations 
has proved difficult (Voigt and Wolf, 1988]. In the i^S 
magnetosphere, Sckopke [1972] studied equilibrium under 
the assumptiM of isotropic plasma pressure in a perturbed 
dipole field. This type of investigation was r 
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tended by Zavnyev and Hasegau* (1989] and Cheng [1992] 
with the more realistic consideration of anisotropy in plasma 
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magnetotad equilibrium has been studied analytically 

[//arm, 1962; Kan . 1973] as well as examined by numeri- 
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^. 1989) : altb ° u ^ 1 thls “ay be merely an artifact of the 
two-dimensional limitation [Kivelson and Spence 1988] 
An alternative approach to obtaining equilibrium con- 
figurations is to use data. For example, tl^gnetic field 

U987 ' 1989] - whicb gilfTvera^ 

magnetic fieW valu« satisfying the static MaxweD’s equa^ 
bons, have been used extensivefy in relating magnetospheric 
regions totte ionospheric level and vicevST^W^ 
t °} h \ 1990; 1991; Elphinstone and 

' 1992) - However ' wi y mod «k of this type’ 

19821 ^S 011 ! 1 °! f ° rCe . balanCe and Southard, 
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(1) What w the plasma pressure distribution in force 
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(2) How is this plasma pressure distribution compared 
with actual measurements? 
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It is quite evident that the availability of plasma pressure 
distribution in force balance with the magnetic forces in an 
empirical field model will further enhance its utilization. 

Spence et al . [1987, 1989] and Kan et al. [1992] have 
inverted the earlier magnetic field models [Tsyganenko and 
Usmanov , 1982; Tsyganenko , 1987] to obtain the equilib- 
rium plasma pressure along the tail axis at midnight by 
assuming isotropic or nearly isotropic plasma. The goal of 
this paper is to extend these previous studies in deriving 
empirical magnetic field models for the quiet time mag- 
netosphere, with the associated plasma pressure in force 
equilibrium with the model magnetic field and in good 
agreement with observations. Two features are incorpo- 
rated, namely, the presence of an eastward inner ring 
current and the occurrence of anisotropic plasma pressure. 
These two features are found to be essential in arriving at 
a realistic plasma pressure profile. This study treats the 
simple situation in which the dipole tilt angle is zero, 
making the configuration symmetric and easy to handle 
without significantly altering the results. 

The Modified Magnetic Field Models 

The magnetic field models of Tsyganenko [1987, 1989], 
hereafter referred to as T87L (long version) and T89, have 
a number of features that lend themselves to easy applica- 
tions. The magnetic field in these models is given in an 
analytical form, with most terms bearing physical insights 
on the various current systems contributing to the total 
magnetic field at a given location. The coefficients of these 
terms are obtained through a least-squares fit to actual 
measurements and are derived separately for different 
levels of the Kp index to allow for the variability of the 
magnetic field with geomagnetic activity. Since we con- 
sider here the quiet time magnetosphere, the activity level 
chosen for this study corresponds to the Kp = 0 level. Modi- 
fied versions of these field models form the basis from 
which magnetic field models are constructed here. 

The ring current is dominated by the gradient of plasma 
pressure. At the inner edge of the ring current, the plasma 
pressure decreases rapidly inward, producing an eastward 
flowing ring current in the inner magnetosphere [Lui and 
Hamilton , 1992]. This ring current is featured in neither 
the T87L nor the T89 model. Although this eastward flow- 
ing ring current may not alter the magnetic field to a great 
extent, as will be demonstrated later, its absence tremen- 
dously affects the plasma pressure distribution required to 
maintain equilibrium with the magnetic forces. Without 
such an eastward flowing current, the plasma pressure in 
force equilibrium with the magnetic field will have to in- 
crease monotonically inward, as indicated by the earlier 
work of Spence et aL [1989] and Kan et al. [1992], and is 
contrary to the observed plasma pressure distribution. Even 
though the required volume current density for the east- 
ward ring current is small, the j X B force is significant, 
since the field strength is high (about inversely propor- 
tional to the third power of the geocentric distance) in the 
inner magnetosphere. 

To include the inner ring current in the model, we con- 
struct its vector potential using a procedure similar to the 
one employed in the T89 model to represent the outer ring 
current, i.e., 

Arc = Qrcp^Trc > (1) 


where C 1RC is a coefficient controlling the intensity of the 
inner ring current and 

S|RC = 'Ip* + (ante + tme) 2 . (2) 

free = V** + + "riRc Arc)* • O) 

Arc = 0.5^1+x{* 2 + Z 1RC , )~ l/a j. (4) 

The location of the point in space is given in cylindrical 
coordinates by (p, z ) and in Cartesian coordinates by (x, 

y, z). The physical significance of the nonlinear parameters 
is the same as given for the outer ring current by 
Tsyganenko [1989]. The radial scale length a mc relates to 
the geocentric distance for the inner ring current. The half 
thickness of the inner ring current region is represented 
by Dq. Tbe variation of this thickness with local time, i.e., 
the difference between dayside and nightside, is controlled 
by the function with its scale length L^. We find 
that the flow reversal of the ring current from westward to 
eastward inside L ■ 3.5 can be reproduced reasonably well 
by adopting the following values for the above parameters: 
^IRC = n# F> ^irc ~ D. = 1.8 R £t y\R£ = -0.2723 

(same as the value for the quiet time outer ring current), 
^d Z/jb £ » 1.5 Rg. To keep the total magnetic moment of 
the ring currents the same as the original ring current in 
the T89 model, we have modified the outer ring current 
strength by 

Crrc = Cor c “ Qrc x ( q irc / °orc ) 2 * W) 

where the subscripts RRC and ORC denote, respectively, 
the revised and original outer ring current parameters in 
T89. It turns out that the modification is extremely slight 
(only -0.2%). Tins model with the inner eastward ring 
current and a correspondingly modified outer ring current 
strength is identified hereafter as T89R. 

Since the T89 and T89R field models give rather weak 
magnetic fields in the magnetotail [Stem and Tsyganenko , 
1992; Rostoker and Skone , 1993; Peredo et al 1993], we 
have also explored the T87L field model. To take advan- 
tage of the fact that the T89R model is better than the 
T87L model in the inner magnetospheric region, but the 
reverse is true for the distant tail [ Tsyganenko , 1989; 
Donovan et al. t 1992], we combine the two field models by 
constructing a hybrid field model (designated hereafter as 
TH) with a transition parameter c such that 

^*,TH “ + (1“ e ^z,89R » (6) 

< = 0.5 - 0.5 sin ((2(x - x L ) /(xy - x L ) - l)r / 2) , (7) 

where B, is the x component of the magnetic field along 
the midnight meridian and subscripts 87L and 89R denote 
parameters from the T87Li and T89R models, respectively. 
The parameters and Xy are the lower and upper bounds 
of x for the transition. The transition parameter € is 0 at 
x = Xy and 1 at x = Xy, with a smooth derivative dddx at 
the end points. The boundary points adopted here are x y 


i 


= 5 R E and x L = 25 R E . We note that the region of 

interest in this paper is the equatorial plane along the tail 
axis in which the magnetic field has only the z component 
(for a dipole tilt angle of 0°). The divergence of the mag- 
netic field along the tail axis therefore reduces to dB Idz. 
It is clear then that the linear combination of the mapietic 
fields from the T87L and T89R models as expressed in (6) 
is also divergence free along the tail axis. On the other 
hand, if this hybrid model were to be extended beyond the 
tail axis, there would be a small nonzero divergence if the 
other field components were not modified accordingly. 

Figure la shows the profiles of magnetic field values in 
logarithmic scale from the four magnetic field models, i.e., 
T87L, T89, T89R, and TH, at the equatorial region along 
the midnight meridian. Note that the four magnetic field 
profiles are quite different between x - -8 R E and -30 R E . 
TTie highest magnetic field values come from T87L, fol- 
lowed by TH, T89, and T89R. Although the inclusion of the 
inner ring current only slightly weakens the magnetic field, 
this effect is manifested quite dramatically in the extremely 
low field region of the T89 model, where the field magni- 
tude goes below 1 nT. 

The volume current densities from these four models at 
the equatorial region along tbe midnight meridian are 
examined in Figure 16. The T87L model typically gives the 
lowest current densities among these models. Tailward of 
* “ —10 R E , the T89 and T89R models give almost the 
same values, whereas the TH model shows noticeably lower 
current densities between x « - 10 and -30 R E . The lower 
values from the TH model, in comparison with the T89 
and T89R models, are expected from the influence of the 
T87L model. Earthward of x - - 10 R E , although the cur- 
rent densities from the T89R and TH models are almost 
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Fig. 1. (a) Magnetic field profile* along the midnight meridian at 
the equatorial plane from four magnetic field models. (6) Current 
density profiles along the midnight meridian at the equatorial 
plane from four magnetic field model*. 


identical, those from the T89 and T87L models have differ- 
ent values. Earthward of x - -8 R £ , both the T89 and 
T87L models show a persistent westward current. In con- 
trast, the T89R and TH models show a reversal of west- 
ward ring current to an eastward ring current at -3.5 R 
The presence of this current reversal is modeled to repre- 
sent the reversal noted from the CCE observations [Lui 
and Hamilton, 1992], In comparison, we find the observed 
current densities appear to be higher than the model val- 
ues at the current density peak but smaller than the model 
values at distances further downstream. In addition, the 
observed current densities are more variable, reflecting the 
filamentary nature of current much like the magnetotail 
current [A fcComas et aL, 1986]. 

The observed anisotropy of plasma pressure taken from 
Lui and Hamilton [1992] is shown in Figure 2 together 
with an empirical fit. The anisotropy is seen to be large in 
the inner region and decreases rather systematically at 
distances further downstream. The empirical fit to the 
anisotropy is a fifth -order polynomial function of x _1 , i.e.. 


^-1 = 00 + «i *' 1 + a?* "* + ogx -3 
+ a 4 x' 4 + o 5 x -6 , 


( 8 ) 


with coefficients Oq = -0.410554, a t = -9.94369, a, = 
-86.9877, a, = -504.066, a 4 = -1110.73, and of = 
-847.912. The range of validity for this fit is -2.5 R E > x 
> ~ 15 R e such that the plasma pressure is isotropic at x 
“ “15 R e . That the pressure should be isotropic down- 
stream of x - - 15 R e is indicated by the studies of Stilez 
et at. [1978] and Baumjohann and Paschmann [1989]. 

Given a magnetic field model, it is straightforward to 
compute the magnetic force j X B and determine the pres- 
sure force required to balance it. The radial profile of per- 
pendicular plasma pressure can be obtained by integrating 

^ = [jxB + (P 1 -P,XbV)b] r , ( 9 ) 


where the subscript r indicates the radial component of 
the vector quantity, b is the unit magnetic field vector 
B, j is the volume current density, and P ± and P, are the 
plasma pressure components perpendicular and parallel to 
the observed magnetic field direction, respectively. 

Figure 3a shows the result from this computation for 
the T89R model with the plasma pressure at x = -34 R £ 
taken to be 0.074 nPa. This initial value of plasma pre^ 
sure is obtained by computing the plasma pressure re- 



Fig. 2. Observed pressure anisotropy profile of the nightside mag- 
netosphere and an empirical fit. 
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Fig. 3 (o) Equatorial profiles of the perpendicular plasma ores 

Mtk “field ^r mght rae, 1 dian inverte<J from the T89R mag- 
®*“ “**'! as fy nun ^ isotropic plasma pressure (dashri 
traoe) and adopting the pressure anisotropy shown in Fieure 2 
(solid trace). (6) Equatorial profile of the perpendicular pUsma 

fi°u 8 .“ ldmght meridian inverted from the T89R 
mapietK: field model in comparison with measured values of per- 
pendicular plasma pressure from several satellites, ** 

quired at the neutral sheet (e = 0) to balance the ‘asymp- 

fu ,d ^ f ° r the ^ ,obe re e° D of 

T89R model at that distance. The ‘asymptotic" magnetic 
field value is the x component of the field at r = 20 R 
which in this case is -13.6 nT. It is worth noting that sin£ 
the plasma pressure is shown in logarithmic scale, a larger 
rniUal value for the plasma pressure will not affect the 
plasma pressure in the inner magnetosphere significantly. 
Ine plasma pressure is assumed isotropic at x = - 15 R 
and further downstream. The computed profile of the per* 
pendicular plasma pressure is given by the solid trace For 
comparison, the profile of plasma pressure, assuming pres- 
aire isotropy at all distances, is given by the dashed trace 
It is quite apparent that the assumption of isotropic pres- 
aire results in an underestimate of the plasma pressure at 
the inner magnetosphere by as much as a factor of -3 

, findlng readil y ^ understood by noting that the 
plasma pressure anisotropy inside * - - 15 R E for the quiet 
time magnetosphere is in the sense of P, > P.. With this 
sense of prepare anisotropy, the second term on the right 
of equation (9) adds to the first one, thus requiring a larger 
pressure gradient for force equilibrium, as pointed out 
earlier by Lut and Hamilton [1992] 

P'"”* p™me j, 

S, «• ATS 5, ISEE 

81x1 CCE m Figure 36. The ISEE 2 data are taken from 

fc," 1 '*S 9 |, f0r .?* * ’ £e 

JTt, Haml t0n 119921 the ATS -5 data from 
Deforest and Mcllwmn [1971], and the Explorer 45 data 

from Spence et al. [1989], using measurements of 1- to 


872-keV protons presented by Smith and Hoffman [1973] 
Pressure values from Explorer 45 data indi- 
cated by the eiror bars actually spans the observed pres- 
sure taken during geomagnetically quiet to disturbed con- 
ditions, and thus the lower limit should be used for com- 

tHe Calculation As can be seen from 
Figure 36, the model pressure profile is in general above 

^ n e n °^. rv ! d , pressure P rofU «. especially in the outer re- 
f_ on - findln g » understandable and is to be expected 
observed pressure i, w „e^r5 SK2 
nght at the equatorial plane and thus is expected to 
underrepresent the actual observed pressure at the equa- 

T InS ', de 4 A P ressure determined from 
“J7 r «« noticeably lower than that from CCE values 
and the model calculation. This difference may be due to 
the upper energy threshold being only 872 keV for Ex- 
fT?? J 5 roeasurements, whereas that of CCE extends to 
4 MeV. It is possible that a significant contribution to 

f^r* re , 4 R P “">* f">» P^fefes above 

® 87 f ke .Y , ^f 1 We have evaluated this possibility using 

We' fou P "“5? ° f 006 I 1 ** 1 

We foun<l the ions with energy of >872 keV contributing 

typically -50-65% to the plasma pressure at geocentr* 

m,de 4 fr Sin0e the ***** m ^ 

nation of pressure between Explorer 45 and CCE is about 
a factor of 6, we conclude that the different energy pass- 
bands can account for a significant portion, but notthe 
entirety, of the observed pressure difference 

th yZ^J e 5^ VCry *”* *«"*““* exists between 
the observed and computed profiles of plasma pressure on 

the basis of the T89R model. Another notable feature of 
the model pressure profile is its relative constancy between 
the downstream distances of 16 and 25 R B . This finding 
indicates the small magnitude of the j x B force in this 
region. Reexamining Figure 1 reveals that the main rea- 
son for the decrease in the j x B force in the T89R model 
w the small value of the magnetic field in that region 

stoSS "“j-* 1 ® V ° ,Ume current density is still sub- 
stantial at those distances (Figure 16). The more abrupt 

m the r ®d“d Profile of plasma pressure occurs in 
the downstream distances between x « -12 R K and x m 
- 15 R s , as noted earlier by Spence et al. [1989], The j X 
B force * therefore relatively large at the transition region 
between the dipolelike and taillike field configure tio^Tlt 
is interesting to note that this region may correspond to 

^ P ° ml thC Uil * Where the tangential stress 
acted on the magnetotail by the solar wind at the tail 

the?!!? T T 1 ** the ®tt r active force between 

the Earth s dipole and the tail current system ISiscoe and 
Cummings, 1969]. 

.. A jj mil “f f lasma pressure profile can be obtained from 
the TH model, as shown in Figures 4a and 46. The plasma 
pressure based on the formula from Spence and Kivelson 

P(nPa) = 89e- 4)59 W + 8.9|xr 163 (10) 

t the downstream distance x 

in Earth radu. The assumption of isotropic pressure 
gives an underestimate of plasma pressure by as much as 
a factor of -3 in the inner magnetosphere as before. The 
most significant difference between the two anisotropic 
pressure profiles in Figures 3 and 4 lies in the downstreL 
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Fig. 4 . (a) Similar to Figure 3o except that the perpendicular 
plaama Pressure* inverted from the TH magneuTneld 
instead. (6) Similar to Figure 36 except that the perpendicular 
plaem. pressure is inverted from the TH magneto field mdW 


region of x- -10 to -20 R g . Whereas the T89R profile 
shows a relatively constant pressure in this region the 
other profile shows a continuous pressure increase toward 
the Earth. This continuous rise, however, only contributes 
rather insignificantly to the plasma pressure earthward of 
x- -8 R s . The continuous change in plasma pressure re- 
flects a more gradual transition between the dipolelike and 
taillike field configurations as well as a more * -dependent 
pressure for the near-Earth tail region in this model than 
that in the T89R model. The empirical formula by Spe nce 
and Kivelson (1993) also compares reasonably well with 
the model profile. The formula gives a good fit to the plasma 
pressure within the range of x - -2.5 to -15 R E but over- 
estimates it inside this region and underestimates it fur- 
ther tailward. 

Summary and Discussion 

We have investigated two empirical models for the speci- 
fication of magnetic field and the plasma pressure for the 
quiet time mghtside magnetosphere. The derived magnetic 
field models are extensions of the previous magnetic field 
models of Tsyganenko [1987, 1989] by the addition of an 
inner eastward ring current. Furthermore, in constructing 
the equilibrium plasma pressure along the tail axis in the 
midnight meridian, we have taken into account the pres- 
ence of plasma pressure anisotropy in the inner magneto- 
sphere earthward of x- -15 R E . An empirical fit to the 
observed pressure anisotropy in the inner magnetosphere 
(between 2.5 and 15 R E ) is also provided. These two new 
features profoundly affect the deduced plasma pressure 
profile. Hie first feature leads to an earthward reduction 
of the plasma pressure in the innermost part of the mag- 


netosphere (geocentric distances of < -3 R E ) as observed 
The second feature leads to an appropriate evaluation of 
the plasma pressure gradient. Without the consideration 
of plasma pressure anisotropy, the deduced plasma pres- 
underestimated by as much as a factor of -3 
The derived perpendicular plasma pressure is found to be 

rrr*'nf! eme jV W ' t1 ' the observed values from ISEE 2 
CCE ATS 5, and Explorer 45. The T89R and TH models 

therefore provide quite realistic representations along the 
tail axis in the midnight meridian of not only the quiet 
tune magnetic field but also the associated equilibrium 
plasma pressure distributions needed to provide the force 
balance with the magnetic forces in these models. The range 
of validity for the T89R and TH models in representing^* 
magnetic field, the volume current density, and the aniso- 
tropic plasma pressure considered here is -2.5 R- > x > 
-36 R e . The present result also supports the earlier find- 
mgs from Spence et at. [1989] that the gradient of plasma 
pressure shows a relatively large change in the transition 
region between dipolelike and taillike field configurations. 

A region called the inner edge of the plaama sheet has 
been introduced in the early studies of the energetic par- 
ticle environment in the magnetosphere Wasyliu na t, 1968- 
Fmn*. 1971]. It is a region characterized by an exponen- 
j. u®crease of electron energy density with decreasing 
radial distance and is generally located at geocentric dis- 
t&nc f! a . 0 * -6 to 8 R s in the nightaide. Many researchers 
implicitly assume that the inner edge of the plasma sheet 
is associated with an earthward decrease of the total plasma 
pressure, and theories on the closure of the large-scale 
2 foW-^gned current system \Jijima and Potemra, 
78J have been built upon this assumption (see, for ex- 
ample, the review by Mauk and Zanetti [1987]). As pointed 
out by Mauk and Zanetti [1987], the observed plasma 
pressure shows a persistent earthward increase rather than 
a decrease in the geocentric distances usually ascribed to 
the inner edge of the plasma sheet. In terms of total plasma 
pressure, there is no edge at those distances. Hie plasma 
commonly identified as constituting the plasma sheet in 
the outer magnetosphere (tailward of the geocentric dis- 
tance of -10 R e ) gradually becomes the hot ring current 
particle population. The absence of a large-scale plasma 
pressure decrease in the so-called inner edge of the plasma 
sheet indicates the necessity to modify the conventional 
mechanism for closing the large-scale region 2 field-aligned 
current. An alternative means for the region 2 field-aligned 
current closure in an earthward increasing plasma pres- 

noQoi°'Tn.‘ tl °v HaS been proposed ^ Ltd and Hamilton 
11992J. They have suggested, based on the work of Sato 

and hjima [1979], that the correct sense of field-aligned 
current will be generated if the VP X vector is inclined to- 
ward the midnight meridian more than the VB vector in 
both the premidnight and the postmidnight sectors, since 
the sense of the field-aligned current is determined by the 
s»gn of the triple product of B, VP X , and VB. Figure 5 
illustrates this configuration graphically. In other words, a 
plasma pressure in the midnight sector slightly enhanced 
over that of the adjacent local time sectors will satisfy the 
requirement 

At least two aspects of this study can be improved in 
the future. One is to extend the calculation to other local 
times besides midnight, and the second is to extend it in 
the z direction to cover the region away from the equate- 
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ABSTRACT 

The complete disappearance of energetic electrons was observed by CRRES in the near geosyn- 
chronous region in 7.5% of the orbits examined. These total flux dropouts were defined by the 
fluxes rapidly dropping to levels below the sensitivity of the MEA energetic electron spectrometer on 
the CRRES satellite. They were separated into those that were only energetic electron dropouts and 
those that were associated with energetic ion and plasma dropouts. Approximately 20% of the events 
showed dropouts of all particle fluxes, and these were usually coincident with large increases in the 
local magnetic intensity and signatures of strong current systems. The energetic particle instruments 
and magnetometer on CRRES provide a detailed picture of the particle and field responses to these 
unusual conditions. Both the local morning and dusk events were associated with strong azimuthal 
(eastward) and radial changes in the magnetic field indicative of a strong current system approaching 
and sometimes crossing the CRRES position at the time of the flux dropouts. The direction of the 
field changes and the details of particle observations arc consistent with CRRES passing through the 
plasma sheet boundary layer and entering the tail lobe for a significant number of the events. 

INTRODUCTION 

Changes in the energetic particle distributions at near geosynchronous altitudes arc often observed in 
association with substorm processes and the large scale boundary motions that often accompany them 
/ 1 /- ^" or exani plc. on the dayside, the compression of the magnetosphere by solar wind pressure 
pulses can move the magnetopause inside the geostationary orbit /2.3/. This is evidenced by the 
dropout" of energetic particle fluxes observed by spacecraft in the dayside magnetosphere /2,3,4/. 
This occurs because the interplanetary energetic particle fluxes are orders of magnitude lower than 
those found inside the magnetosphere. Thus, when the boundary is "pushed" earthwards of a mag- 
netospheric satellite's orbital position, it changes from observing high magnetospheric fluxes to 
observing the low interplanetary flux in a very short time (i.e., the flux drops out). It is relatively 
uncommon to observe such magnetopause crossing events near geosynchronous orbit /2 /. In the 
nightside magnetosphere, the distortion of the geomagnetic field caused by the intensified cross-tail 
currents which arise during the substorm growth phase can cause significant modification in the 
energetic particle angular distributions and intensities. These can become so extreme as to cause 
complete flux dropouts of the energetic particles (see for example refs. /5/ and 161 ). Normally, such 
signatures are observed in the midnight sector, and often the fluxes only decrease drastically but do 
not totally drop out. If magnetospheric satellites enter the tail lobe, which does not contain signifi- 
cant fluxes of energetic panicles under normal conditions, the energetic panicles would “drop out". 
Normally, near-earth satellites, such as those with orbits that do not extend beyond geosynchronous 
altitudes, do not enter the tail lobes. Often the post-dropout recovery of the energetic particle fluxes 
are associated with the expansion onset of a substorm /5/. It is relatively uncommon to observe strong 
energetic panicle flux dropouts in the dawn and dusk sectors /l/. Korth et al. /7/ has recently reported 
the occurrence of energetic particle dropouts in the local morning sector and Moldwin et al. /8 / have 
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observed flux dropouts at nearly all local limes and tail lobe entry by geosynchronous satellites from 
dusk to dawn in the nightside magnetosphere. In this paper we use data similar to that of Korth et al. 
PI but extend the study to include ail local times reached by CRRES (Combined Release and 
Radiation Effects Satellite). By using the local magnetic field and particle data in combination we 
can show the relationship between the particle dropouts and the active current systems at the time of 
such events. 


INSTRUMENTATION 

The CRRES satellite was launched on 25 July 1990 into an elliptical orbit with a geocentric perigee of 
-<5720 km. apogee of -39950 km and an inclination of 18.15°. /9/ The initial local time of apogee 
was -8.7 MLT. and after 1072 orbits (near end of life 14 October 1991) it was at -14 MLT. The 
orbital inclination in combination with its -9.6 hour orbital period and the Earth's tilted magnetic 
dipole allowed CRRES to reach L shells up to and somewhat beyond L=8 on a regular basis. The 
CRRES near apogee altitude coverage has some gaps in local time. These are the local time intervals 
near noon on the dayside (from —09 through 13 MLT) and in the early post-midnight region near 01 
MLT. The post-midnight coverage gap resulted from the suspension of instrument operations during 
the CRRES near-apogee eclipse periods. 

The instruments used for this study are primarily the Medium Electrons A (MEA) Spectrometer /10/, 
the Magnetosphere Ion Composition Sensor (MICS) /ll/. and the CRRES magnetometer /12/. We 
have also examined some EPAS (Electron and Proton Wide-Angle Spectrometer) data /1 3/ and the 
LEPA (Low Energy Plasma Analyzer) data /14.15/ to determine whether a plasma signature was 
present or not. The MEA measures electrons from -153 keV to 1.5 MeV in 17 differential energy 
channels /10/, while the MICS measures ions from -1.2 to 426 keV/charge with ion charge and mass 
composition over mass dependent energy ranges /ll/. The MEA and MICS sensors are mounted 
perpendicular to the CRRES spin axis, which is sun pointed at all times, and cover a relatively wide 
range of particle pitch angles during a spin period (-30 sec). The EPAS sensor covers the energy 
ranges of 21 - 285 keV and 37 - 3200 keV for electrons and ions, respectively /13/. EPAS has mul- 
tiple fields of view (FOV) and, by combining the data from the multiple FOV, basically obtains a 
complete pitch angle distribution every spin period. The LEPA also has multiple FOV and provides 
complete pitch angle coverage of electrons and ions with energies of 0.12 to 28 keV/q. In the LEPA 
summary data referenced here /1 5/, only the precipitating and perpendicular plasma electron and ion 
fluxes were available. 

OBSERVATIONS 

For this study we required that the one-minute spin-averaged MEA fluxes, at all energies > 153 keV, 
drop to background levels. This is a rather stringent requirement and basically deletes most of the 
growth-phase-only signatures /5/ from the study. The MEA criteria were selected because of the easy 
access to the data and the large geometric factor and good sensitivity of this instrument (see Table I). 
The total MEA flux dropout to background levels was chosen to eliminate the more usual flux 
decreases that occur commonly in the magnetosphere near geosynchronous altitudes. A total of 1 17 
flux dropout events met this requirement between CRRES orbits 76 and 1065 (from 26 August 1990 
through 11 October 1991). There were as many as five separate flux dropouts observed in a single 
orbit. In total, the 1 17 dropouts were distributed over 70 orbits, as shown in Table 2. 

The 117 events were further reduced to those which showed a dropout of the energetic ions 
(>20 keV). Only 53 events satisfied this criterion. Finally, we selected only those events which 
showed a dropout of the plasma ions and plasma electrons, respectively, as shown in Table 2. From 
the initial selection of 1 17 events, only 20 showed a dropout of all particle fluxes. 

Flux dropouts were observed to occur at all local times on the nightside magnetosphere. (We 
specifically eliminated obvious cases of magnetopause crossings from this study.) The most unex- 
pected regions for them to occur are the local morning and local evening. In these regions the mag- 
netic field topology is usually fairly dipole-like for L shells below 7. Thus, we expect that the 
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distorted magnetic configurations that lead to flux dropouts should be rarely observed there by 
CRRES, even at the moderate latitudes it attains. Also, one would naively expect that the local 
morning and local evening occurrences should show similar features. We present below one case 
each of local evening and local morning flux dropouts as examples of the phenomena observed by 
CRRES before we discuss the statistics of the dropouts. It should be noted that CRRES was in the 
northern hemisphere near local morning apogee and in the southern hemisphere near local evening 
apogee. 6 

TABLE 1 . Instrument Parameters 


MEA 


Geometric Factor 
2.14-5.88 cm 2 srkeV 


MICS 


[ref. /10/] 

4.4 x 10 cm 2 sr keV 
[ref. /!!/) 


LEPA 


Minimum M exsurib lt Flux 
~ 3 x I0‘ 2 Eiectrons/(cm 2 sec sr keV) 


- 23 Ions/(cm 2 sec sr keV) 

- 5 x 10* p*rticles/(cm 2 sr) »i 0. 1 2 keV 
(estimated from summary spectrograms] 


TABLE 2- Particle Flux Dropout Statistics 


1 

Number of 

Kp Value 

IDstl Vxlue 

Dropouts 

Orbits 

<3 

3-6 

> 6 

< 30 

30 - 100 

> 100 

Energetic Electrons 

117 

70 

17 

65 

35 

34 

46 

37 

Energetic Ions 

53 

34 

6 

27 

20 

13 

21 

19 

Plasma Ions 

33 

20 

0 

20 

13 

5 

15 

13 

Plasma Electrons 

23 

14 

0 

17 

6 

4 

10 

9 

AH Particles 

20 

VJ _ e _ 

14 

0 

15 

5 

4 

10 

6 


Number of orbits examined = 944 


Local Evening Flux Dropout 

An example of an event which showed the dropout of all particle fluxes on the local evening side of 
the magnetosphere is presented in Figures 1 and 2. Figure 1 shows the energetic electron data from 
MEA for this event, which occurred during CRRES oibit 767 on 5 June 1991. The flux from several 
different electron energies are shown on the plot, with the highest trace being that of the - 153 keV 
electrons. The decreasing flux levels correspond to increasing energy. The bold trace represents the 
response of a "background" monitor (BKG) within the MEA. The -153 keV electron fluxes were 
required to drop to or below the BKG level before they were identified as a total flux dropout There 
were three such periods of total electron flux dropout identified in Figure 1: one near 1500 UT at 
L=7.0, 17.8 MLT and -21.8° magnetic latitude (MLAT). the second near 1600 UT at L=7.4, 18.5 
MLT and MLAT = -21.6°. and the third near 1850 UT at L=6.0, 20.7 MLT and MLAT = -18.5°. 


The lower panels of Figure 2 show energy-time spectrograms of several channels of MICS data for 
the same orbit. DCR is a measure of the total ion flux independent of species and thus is similar to 
the measurement made by a total ion plasma instrument such as the LEPA. only at somewhat higher 
energies (- 7 - 400 keV/q versus 0.12-28 keV/q). The MICS Alpha" channel represents the total 
flux of all ions with mass > 4 AMU. The other channels measure the fluxes of the H + , He + , He ++ 
and 0+ above a mass dependent energy threshold /1 1/. There are periods in the MICS data where 
essentially all ions dropped out, especially during the 1600 UT event Note that each 0+ spectrum 
has been averaged over 8 minutes, which means there were a few counts in each interval that straddles 
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Fig. 1. Energetic electron fluxes from the CRRES MEA sensor during orbit 767 on 5 June 1991. The 
topmost trace is the flux of -153 keV electrons and the succeeding traces are for successively higher 
energies /9/. The bold line represents the response to a background monitor (BKG). The -153 keV 
electron fluxes were required to be at or below BKG in order to qualify as a flux dropout for this study. 
For this period, the energetic electron flux dropouts occurred at 54100, 57600 and 67800 sec (or 1500 1600 
and 1850 UT). 
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the flux dropout’s edges. It should be noted that the He + flux dropped well before the rest of the 
ions and that it recovered later than the other ions too. This relationship between the two charge 
states of Helium was often observed in the flux dropout events. This would be expected if He ++ is 
presumed to be the dominant He ion in the outer magnetosphere, as is expected for a solar wind 
source in which He + is generated via charge exchange from He ++ . The He + would be most prevalent 
on the lower L shells deep inside the plasma sheet. If the flux dropout were the result of a 
reconfiguration of the magnetosphere that caused the CRRES satellite to be on field lines that thread 
the distant tail near the interface between the plasma sheet and tail lobe field lines, then the ion com- 
position just prior to and just after the flux dropout should reflect that of the distant plasma sheet, as 
was apparently the case for this event. 

The upper three panels of Figure 2 show the magnetometer measurements for this event The data 
are presented as the differences between the Tsyganenko 1987 field model /16/ for Kp = 3 and the 
observed values. The direction ( b ) was taken as the Tsyganenko field direction. The reference 
coordinate directions b, r and w are in the directions parallel to the Tsyganenko field, earthward in 
the plane containing band the radius vector from the center of the earth to CRRES and westwards 
respectively, as shown in Figure 3. 

5Bb, 9Br and 9Bw are the differences (residuals) between the measured and the Tsyganenko model 
field. As can be seen by the magnitudes of these components, especially 9Br. the magnetic field 
intensity at CRRES was stronger than normal and very "tail-like”. There was also a significant but 
varying, westward component. The large 9B r and 9B W indicate that strong currents were flowing in 
the neighborhood of CRRES. Specifically, near the onsets of the particle dropouts, the field became 
more tail-like and rotated from a westward to an eastward direction. At each flux recovery the field 
became less tail-like (although still much more tail-like than normal) and rotated back towards the 
meridian plane (where 9Bw = 0). In particular, a strong dipolarization occurred in conjunction with 

?^>n 0 . n r^ UX u reC0Very " ear 15,5 07 A second (momentary) and third dipolarization occurred near 
620 UT, where there was first a short ion flux recovery with subsequent dropout followed by a per- 
manent ion flux recovery near 1628 UT (see DCR in Figure 2 near 15200-1530 UT). These field 
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Fig. 2. CRRES magnetometer and MICS data for 11-21 UT on 5 June 1991 Fig. 5. Like Figure 2 except 09-19 UT on 7 September 1990 (orbit 106). Top 

(orbit 767). Top three panels show the deviations of the three components of three panels show the deviations of the three components of the magnetic field 

the magnetic field from the Tsyganenko model /\ 6 /. The bottom six panels from the Tsyganenko model /\ 6 /. The bottom five panels show energy-time 

show energy-time spectrograms of different ion species measured by the spectrograms of different ion species measured by the MICS instrument (see 

MICS instrument (see text). text). 
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signatures are consistent with the CRRES satellite having been approached and crossed by a Region 1 
current system as the particle fluxes dropped out and then recovered. 

While there was a total energetic electron flux dropout near 1850 UT, there was not a comparable 
dropout in the energetic ion fluxes. The intensification of the ion fluxes near 1700 UT (61200 sec) 
occurred simultaneously with a weak energetic electron recovery. The ion fluxes returned to lower 
levels by 1730 UT, which is just prior to the electron flux recovery that peaks near 1752 UT. 

The plasma data (not displayed here) showed a total plasma electron dropout in conjunction with 
both of the energetic ion dropouts in Figure 2. The plasma ions also dropped out but only just prior 
to the energetic particle flux recoveries near 1515 UT and 1625-1628 UT. Both the plasma electrons 
and ions showed the same momentary flux recovery, subsequent short dropout and final recovery 
during this latter period, as was observed in the energetic ions (ref. Figure 2). Like the energetic 
ions, the plasma data did not show a flux dropout near 1850 UT and indicated that CRRES remained 
in a hot plasma sheet plasma. The only other significant feature in the plasma data was the occur- 
rence of low energy (< 1 keV) field-aligned electron fluxes just prior to and immediately after the 
plasma and energetic particle dropouts discussed above. 

Local Mo rning Flux Dropout 

Figure 4 shows an example of energetic electron total flux dropouts observed by CRRES near local 
morning on 7 September 1990. There were five flux dropout intervals (marked by shaded vertical 
bars in Figure 4) which started near 1300. 1308, 1320, 1332 and 1409 UT. We have plotted the 
CRRES magnetic difference data (as discussed above) on the same plot for comparison. Note that the 
dBr was positive when the magnetic field had tail-like conditions because CRRES was in the northern 
hemisphere, whereas 5Br was negative during the local evening tail-like conditions discussed above. 


N 

<» 


Fig. 3. Coordinate system used for detrended magnetic 
field data. £ is parallel to the Tsyganenko field vector, r 
is earthward in the plane defined by l> and the radius 
vector from the center of the earth to CRRES and w is 
westward. 



Fig. 4. Energetic electrons and detrended magnetic 
field data for 7 September 1990 (CRRES orbit 
106). The flux dropout onsets are 1300, 1308. 
1320. 1332 and 1409 UT. 
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In both cases, the dBw turned from westward towards zero or eastward at the flux dropout onsets and 

!^^STJ rd - at ‘Si h UX reCOvery - Likc ^ local evcnin * observations, the field is compressed 
(positive 5Bb) dunng the dropouts. Also, the field was tail-like at the dropout onset and dipolarized 

* reC °y!7 Thcse * ,gnan,res> ,ike of *e local evening dropouts, are consistent with 
CRRES having been approached and crossed by a region 1 current system at the flux dropout and at 
recovery. This would be consistent with CRRES having entered the tail lobe on the local morning 
side of the magnetosphere. 6 


Figure 5 shows the MICS energetic ion data and, again, the magnetic field difference plots for this 
event. In this case, there was a one-to-one correspondence between the energetic electron flux 
dropouts and the energetic ion flux dropouts. The primary difference between Figure 5 and Figure 2 
was the relative absence of 0 + ions throughout the magnetosphere on 7 September . Other than this 
lack of 0 and He ions on 7 September, flux dropout features were remarkably identical for these 
two events that were observed on opposite sides of the magnetosphere. 


Examination of the plasma data (not shown here) showed that there was a one-to-one correspondence 
between the energetic panicle flux dropouts and the plasma dropouts for the 7 September events 
The narrow dropouts at 1300 and 1308 UT were barely resolved in the summary plasma electron 
p ots and somewhat less clear in the ion plots. Like the local evening case, intense low energy field- 
igned electron fluxes were observed just prior to and immediately after the flux dropouts. In 
general, such field-aligned "soft” electron fluxes were observed near the majority of the energetic 
ion flux dropouts and at nearly all of the plasma dropouts. 


Flux Dropout Statistics 


As mentioned above, we examined the energetic ion, magnetic field and plasma data that were avail- 
a e or all 117 energetic flux dropout events. Table 2 summarizes the breakdown into the different 
categories according to whether the energetic ions and plasma showed flux dropouts and some crude 
ranges of K P and D S x levels for the corresponding events. In general, the K P values observed at the 
time of the flux dropouts were moderate to high. This is presented graphically in Figure 6. where the 
distributions of both energetic electron and energetic ion flux dropouts in K P and Dst are 
shown. The dominance of moderate to high K P values indicates the flux dropouts occur during 
disturbed periods, as would be expected for the strong currents observed in association with the 
dropouts (e.g.. see Figures 2 and 5). But, there is no strong dependence on Dgj The majority of the 
flux dropouts occurred during relatively low or modest D ST This is somewhat surprising. It 
indicates that magnetic storm conditions are not necessary for the flux dropouts to occur. One might 
have expected that flux dropouts would have been more prevalent during the magnetically disturbed 
conditions generally associated with larger ring current enhancements. 


Figures 7 and 8 show the distribution of the CRRES flux dropouts in magnetic local time and L and 
in magnetic latitude and L. respectively. Note that the paucity of events in two different local time 
regions, one post-midnight (0 - 2 MLT) and one near noon (9-14 MLT), partially reflects the lack of 
coverage in these regions. [CRRES apogee never reached the near noon sector and the instruments 
were turned off in the post-midnight sector during the eclipse season.) There was only one obvious 
magnetopause crossing observed by CRRES and this was eliminated from the data set. The different 
panels in Figures 7 and 8 correspond to the distribution of energetic electron dropouts (a), the ener- 
getic ion dropouts (b) and the plasma ion dropouts (c). respectively. Note that when one required 
that the plasma ions disappear to have a total panicle dropout, there remained only two events that 
met this criterion on the morning side but a significant number on the evening side. All flux dropout 
events were observed at magnetic latitudes 10° or more above and below the magnetic equator. The 
preponderance of the dropouts occurred for L > 6.6 and all but two occurred for L > 5.5. This 
indicates that to have a good probability of observing panicle dropouts requires a satellite to be at or 
above geosynchronous altitude and/or at moderate to high magnetic latitudes. 


SRFs n irhL t i 1C rJ^r ntS inFigu ": 8 t0 n J ove t0 hi 8 her latitude with larger L shell is a result of the 
CRRES orbital configuration. As noted in the introduction. CRRES only reached L < 6.3 at the 
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Fig. 6. Energetic ion and electron dropout occurrence frequency versus Kp and D<jj. 

magnetic equator, but because of its moderate 18° inclination it could reach higher L shells when 
apogee was off the magnetic equator. The -10 hour orbital period allowed the CRRES apogee to 
pass through a significant range of magnetic latitudes on a weekly basis. The morning side data were 
taken in the late summer and fall of 1990 while the midnight to dusk data were taken in the spring 
through summer of 1991. The combination of Figures 7 and 8 show that the morning side events all 
occurred in the northern hemisphere and the midnight through dusk events occurred in the southern 
hemisphere. 

Finally, we looked at the changes in the local magnetic field at the times of the total plasma dropout 
onsets to see if they were all consistent with the picture, described above, of the CRRES satellite 
crossing (or being crossed by) a Region 1 current system, entering the tail lobes and then returning at 
flux recovery. We grouped the events into morning side and dusk side local time regimes, as shown 
in Table 3. We examined in detail the pre- and post-dropout field trends and the changes that 
occurred at flux dropout onset and recovery. We categorized the changes, as shown in Table 3, 
according to whether the field was tail-like (had a significant radial component) or was becoming tail- 
like at onset and whether the azimuthal component changed and in which direction. In all cases, the 
field was either strongly tail-like or became tail-like at flux dropout onseL In the majority of cases 
the azimuthal component turned eastward at onset Since all the morning side events occurred in the 
northern hemisphere and the dusk side events were in the southern hemisphere, an eastward rotation 
of the field would be consistent with crossing a Region 1 current system in the respective local time 

S ° 11 * S Clear fr ° m Table 3 1,131 *** ma j° ri, y of total plasma dropout events occurred as the 
CRRES satellite crossed (or was crossed by) the boundary between the plasma sheet and tail lobe as 
evidenced by the passage of the Region 1 current system along with the strongly tail-like field 
geometry and the loss of all particle fluxes. 







Particle Flux Dropout Events 

TABLE 3- Field Changes During Total Plasma Dropout Onsets 
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Energetic Electron Total Flux Dropouts 
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Fig. 7. Distribution of flux dropouts in magnetic local time and L shell: (a) for engergetic electrons, (b) for 
energetic ions, and (c) for plasma ions. 
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Fig. 8. Distribution of flux dropouts in magnetic latitude and L shell: (a) for energetic 
electrons, (b) for energetic ions, and (c) for plasma ions. 
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SUMMARY AND CONCLUSIONS 

The local morning and local evening flux dropout events described here and the total plasma 
dropouts represented in Table 3 clearly show that CRRES entered the magnetotail lobe during these 
flux dropouts. In general, as seen in the two events shown, the recovery from the flux dropouts was 
n& associated with a signature of substorm expansion onset at CRRES. Given that most of the events 
were observed away from the midnight sector, it is possible that the substorm signature was not 
observable by CRRES. To determine whether there was a substorm associated with these events 
requires examining ground magnetometer and other satellite data, which is beyond the scope of the 
present work, but will be done as part of a continuing effort on this subject (see Korth et al. /If). 

We emphasize that care was taken to eliminate magnetopause crossings (only one found) from this 
data set by examining the plasma data for signatures of magnetosheath plasma. In all cases, very 
significant currents were flowing near the CRRES satellite as evidenced by the strong changes in the 
magnetic field topology from average conditions. The field was generally tail-like, indicating that the 
satellite was near the inner edge of the cross-tail current system. The actual current system configu- 
ration that can give rise to such strongly tail-like conditions as were observed near 5.5-6.S MLT (see 
Figure 5) and 17.8-19.0 MLT (see Figure 2) is not clear, and its details require more observational 
constraints than are provided by a single satellite measurement 

It should be noted that the majority of the energetic particle flux dropouts included in this data set 
had rapid onsets like those shown in Figures 1, 2, 4 and 5. The criterion that the lowest energy 
channel from the CRRES MEA sensor must drop to background levels removed most of the events 
that had slowly decreasing growth-phase-like flux decrease signatures /5/. While this reduced the total 
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number of electron flux dropout candidates, it provided a list of those most likely to also show 
energetic ion and plasma dropouts. The high sensitivity of both the energetic electron and ion 
sensors (ME A and MICS. see Table I) allowed us to set very stringent minimum flux levels for our 
dropout identifications. As can been seen in Figures 1. 2. 4 and 5, these thresholds were orders of 
magnitude below typical magnetospheric levels and comparable to the cosmic ray background 
response of the sensors. It should also be noted that all particle data in this study were averaged over 
one or more satellite rotations so that details of the panicle angular distributions were not visible. 

The difference between the number of energetic electron and energetic ion flux dropouts is most 
likely a result of the .ons greater gyro-radii. That is. the ions are observable within an ion gyro- 
radius of the plasma-sheet/tail-lobe boundary on the tail lobe side of the boundary, which is much 
deeper in the tail lobe than an energetic electron can reach. In a follow-on study we will examine the 
energetic ions and electrons in greater detail and use the available angular distribution data /1 3/ to 
sound the boundary in the neighborhood of the flux dropout intervals. This should help to more 
clearly delineate the differences between the energetic electrons and ions during the flux dropout 
events. It is most likely that, in some cases, the satellite approached to within a particle gyro-radius of 
the boundary but did not pass through it. 

The local time and seasonal results described above are in apparent disagreement with the recent 
results of Moldwm et al. for events observed at geosynchronous orbit, wherein they found a prefer- 
ence for the southern hemisphere dropouts to occur during the winter season and the northern hemi- 
sphere dropouts to occur during the summer season /8/. In our case, the season did not appear to 

BU k- ^ Mo . ld J win et ^ results were from data taken ncar geographic equator, whereas the 
<~ddcc °K bU T?* 3 mUCh WldCr range 0f 8 eo 8 ra pbic and geomagnetic latitudes. Therefore, 
CRRES observed the dropouts only when it was more than 10° away from the geomagnetic equator 
We d.d not note the geographic position of the dropouts during this study. This will be done during 
the expanded study noted above. 6 


At this point we do not know whether the field-aligned plasma electrons, observed at the edges of the 
ux dropouts, are moving towards or away from the earth since the actual direction is not indicated in 
the summaiy data. This will be examined in greater detail in a later paper. It is important because it 
will more clearly identify the plasma regime that CRRES was immersed in just prior to the flux 
dropouts and at the recovery. Given the strong field-aligned currents detected by the magnetometer 
in these regions, it is possible these intense electrons may be the current carrying particles. 

It IS clear from the above discussion that much remains to be done on the subject of flux dropouts 
Some of the questions that arise are: How can the field be so tail-like at dusk and dawn on such low 
s e s (< 7)? Why are the electron dropouts more prevalent than ion dropouts? What are the 
motions of the tail-lobe boundary relative to CRRES at the times of the dropouts? What is the 
configuration of the tail-lobe boundary and the corresponding current systems during such events? 
Were the solar wind conditions extreme (i.e., was the magnetosphere compressed)? Were substorm 
onsets associated with many of the events? Were such flux dropout events observable at other 
spacecraft, such as those in geosynchronous orbit. The answer to these and other questions will be 
the focus of our future work on this subject. 
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GREAT DEBATES 


IN SPACE PHYSICS 

The What, Where, When, and 
Why of Magnetospheric 
Substorm Triggers 


Harlan E. Spence 


Over the past three decades, terrestrial 
magnetospheric physics has had a unifying 
and hotly debated focus: the magneto- 
spheric substorm. A magnetospheric sub- 
storm is a three-phase phenomenon 
[McPherron, 1979] in which energy is first ex- 
tracted from the solar wind flow, transported, 
and stored within the Earth’s magnetospheric 
magnetic fields (growth phase). The stored 
magnetic energy is then converted and re- 
leased explosively within the magnetosphere 
and auroral ionosphere (expansion phase). 
Then the magnetosphere and auroral iono- 
sphere relax, entering a quiescent state (re- 
covery phase). 

Understanding of this global process has 
unified the magnetospheric and auroral com- 
munities by providing the "big picture" of 
magnetospheric dynamics as a backdrop 
against which unrelated areas of more fo- 
cused research may be put in context. The 
topic debated in the accompanying articles 
centers on one aspect of magnetospheric sub- 
storms, namely what triggers the expansion 
phase. The companion papers present two 
perspectives on substorm triggers: one advo- 
cates a trigger that is driven by an external 
change in the interplanetary magnetic field 
(IMF), while the other argues that external 
triggers are just one of many mechanisms 
that lead to the expansion phase, including 
internal instabilities. 

What? 

To appreciate the importance of the sub- 
storm "trigger," it is imperative to establish 
what a magnetospheric substorm is, at least 
phenomenologically. A consensus definition 
proposed many years ago stated that "A mag- 
netospheric substorm is a transient process 
initiated on the night side of the Earth in 
which a significant amount of energy derived 
from the solar wind-magnetosphere interac- 
tion is deposited in the auroral ionosphere 
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and magnetosphere" [Rostoker et al , 1980], 
The transient substorm effects are global, 
producing dramatic changes throughout 
geospace. These include major topological 
changes in magnetotail geometry' as growth 
phase magnetic field stress is released, ther- 
malization and bulk acceleration of magneto- 
spheric particle populations both in the inner 
magnetosphere and the magnetotail as mag- 
netic energy is converted to particle energy, 
and dramatic intensifications and motions of 
auroral emissions as electrical currents are 
driven in the auroral ionosphere. 

The agreed upon definition above does 
not provide any physical explanation for 
what causes a magnetospheric substorm. A 
complete description requires a detailed ex- 
planation of how the solar wind flow cou- 
ples energy into and through the 
magnetospheric system. Discussions of this 
subject have been the root of most contempo- 
rary debates in substorm physics. To keep 
study focused on what triggers the substorm, 

1 avoid endorsing any one substorm model. 

Two articles in this issue of Eos debate the 
following question: Are substorms triggered 
externally or internally? Lawrence Lyons, a 
scientist from The Aerospace Corporation, ar- 
gues in favor of an external trigger while An- 
thony Lui, a scientist from the Johns Hopkins 
University Applied Physics Laboratory, ar- 
gues for an internal trigger. Both authors pre- 
sent evidence to support their opposing 
views. To put their arguments in perspective, 
it is useful to first review the related body of 
work that precedes them. 

If substorm onsets are triggered exter- 
nally, phenomenological relationships be- 
tween onsets and external boundary 
parameters should exist. Studies that investi- 
gated these factors were a thrust of early sub- 
storm research. Magnetospheric substorm 
strength (measured by the intensity of cur- 
rents flowing in the ionosphere) was found to 
correlate strongly with the magnitude and 
sign of the north-south component of the in- 
terplanetary magnetic field (IMF). In most 
models, solar wind flow energy is imparted to 
the magnetosphere by way of dayside mag- 
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netic field merging between the southward 
component of the IMF carried past the Earth 
by the solar wind, and the northward di- 
rected geomagnetic field near the subsolar 
magnetopause boundary. In this picture, the 
coupled magnetosphere-ionosphere system 
is energized by its interaction with the solar 
wind/IMF. It is most highly energized when 
the IMF is most strongly southward. This 
strong control exerted by the IMF on sub- 
storm strength gave rise to the earlier related 
debate: Are substorm currents driven directly 
by the solar wind/IMF interaction (i.e., exter- 
nally driven) or does the magnetosphere go 
through an energy loading-unloading proc- 
ess (which consists of external and possibly 
internal aspects)? As is the case in most 
"black-versus-white" arguments, subsequent 
studies demonstrated that the answer is 
"gray." Both driven and loading-unloading 
processes contribute to the magnitude of 
magnetospheric substorm currents. 

Study next focused on whether the IMF 
controls not only the magnitude of sub- 
storms, but also directly triggers substorm on- 
set. That brings us back to our main 
discussion. In June 1994, 1 cochaired a ses- 
sion devoted to substorm triggers with Terry 
Onsagerof the NOAA/Space Environment 
Center at a National Science Foundation 
(NSF) Geospace Environment Modeling 
(GEM) Workshop in Snowmass, Colo. At this 
session, "Tail Geometry: IMF Control and Trig- 
gers," we asked Robert McPherron of UCLA 
to review studies from the 1970s and early 
1980s that sought relationships between sub- 
storm expansion phase onsets and variations 
in the direction of the IMF (particularly north- 
ward excursions) or solar wind variations 
(for example, dynamic pressure changes dis- 
cussed in Lui’s article). McPherron reported 
a 45% likelihood that northward turnings of 
the IMF occur at the time of sharp substorm 
onsets, but that better determination of on- 
set time and IMF data are needed. Workshop 
participants agreed that this was a topic that 
could be better quantified with both existing 
data sets and those to be obtained from up- 
coming solar-terrestrial missions. At this 
year’s NSF/GEM summer workshop, Larry 
Lyons reported on a reinvestigation of the 
IMF-substorm connection. He made the 
strong assertion that "...most, and perhaps 
all, substorm expansion phases are triggered 
directly by IMF changes." This conclusion is 
the basis of his new view of substorms 
[Lyons, 1995]. That some or even many sub- 
storms are triggered by northward turnings is 
not controversial, but Lyons’ statement that 
most or perhaps all fit this category, drew the 
attention of many people. Lyons recategor- 
izes geomagnetic activity that lacks an IMF-in- 
duced trigger and that has previously been 
identified as a substorm to a different class of 
phenomenon. This reclassification yields his 
strong conclusion. That redefinition is a mat- 



ter of ongoing debate and emphasizes the dif- 
ficulty with semantics. 

Tony Lui, who explains the otherside of 
this debate, is an advocate of substorm trig- 
ger models based on plasma instabilities 
within the magnetosphere. In such models, a 
local instability occurs within the magneto- 
sphere, perhaps as the result of global bound- 
ary conditions established by the solar 
wind/IMF. The internal instability however 
may occur without any variations in the IMF 
directly. Luis article explores pheno- 
menological evidence that supports the inter- 
nal trigger hypothesis. Rather than searching 
for explicit signatures of an internal trigger, 
which may vary according to the chosen 
physical model, he seeks evidence for the ab- 
sence of a consistent IMF variation associ- 
ated with substorm onsets. Lui’s review of 
the literature suggests that substorms can be 
triggered by a variety of external stimuli, but 
not exclusively by northward turnings of the 
IMF. This conclusion appears to be more 
widely accepted. 

External and internal substorm trigger 
sources are entirely different scenarios that, in 
an ideal world, should be readily discernible 
and could be used to test the strengths of Lyons’ 
and Lui’s conclusions. Why, then, has this de- 
bate continued? The answer may in part be 
merely semantics. That aside, as we have 
learned in other fields of geophysics where na- 
ture is an uncontrolled "laboratory," for exam- 
ple, in earthquake prediction, that identifying a 
trigger mechanism is not always as straightfor- 
ward as it appears. A shortage of possible physi- 
cal theories is not the problem. Identifying 
trigger mechanisms is complicated by the avail- 
able data. Imperfect observations are some- 
times taken at the right times or places, or the 
correct measurements are obtained but at the 
wrong times or places. If we cannot determine 
trigger location and time well enough, any corre- 
lation with associated solar wind/IMF variations, 
which have their own measurement uncertain- 
ties, becomes problematic and leads to unclear 
answers. Uncertainties, then, in the "where" and 
"when" of the substorm trigger have limited our 
ability to reveal the "why." 

Where? 

One important aspect of the physics of the 
substorm trigger is the location where it oper- 
ates. While the growth, expansion, and recov- 
ery phases all exhibit global behavior during 
substorms, it is not as clear where the expan- 
sion phase onset begins, initiating the global 
reconfiguration. Substorm phenomenology 
has very different "faces" depending upon 


where a particular measurement is made 
(that is, ground-based, ionosphere, geosta- 
tionary orbit, or magnetotail) and what quan- 
tity is being measured (that is, magnetic or 
electric fields, charged particles, auroral lu- 
minosity, etc.). These factors influence the in- 
ferred onset location. Consequently, location 
estimates span a broad range, from near-geo- 
stationary orbit to the distant magnetotail to 
the auroral ionosphere. The size of the mag- 
netosphere further confounds a simple reso- 
lution to this question. It is the largest of the 
geospheres with a volume of ~10 3 ‘cm 3 
Spacecraft sample only a minute portion of 
this volume at any one time and, owing to the 
magnetosphere’s dynamism, often cannot 
uniquely distinguish spatial from temporal 
variations. In spite of the enormity of the 
volume, the sparseness of data, and time de- 
pendence, heroic efforts have been made in 
our understanding of the magnetosphere 
and substorms. This is truly a triumph of 
space physics. Critical stages have included: 
establishing global auroral substorm phe- 
nomenology based on syntheses of local, 
ground-based observations [ Akasofu , 1964]; 
quantifying magnetospheric boundaries and 
the different regions they delineate from ex- 
ploratory spacecraft missions; and investigat- 
ing the physics of individual magnetospheric 
regions using focused spacecraft missions. 
New capabilities provided through the coor- 
dinated, multipoint measurements of the In- 
ternational SolarTerrestrial Physics Program 
(1STP) will allow us to advance our growing 
knowledge of substorm onset location. 


When? 

The question of where the trigger occurs 
is intertwined with the question of when sig- 
natures at various locations occur relative to 
one another. The magnetosphere-iono- 
sphere system is electromagnetically cou- 
pled, so despite the great volume, substorm 
dynamics are evident across vast regions on 
relatively short timescales (seconds to min- 
utes). As noted above, the differentiation of 
spatial and temporal variation is an inher- 
ently very difficult task, and especially so for 
substorm studies where strong spatial gradi- 
ents and temporal evolution are both impor- 
tant. The question of when the various key 
phenomena occur at different locations and 
how they connect in a physical sense is 
therefore unresolved. High time resolution 
measurements throughout the critical active 
regions are still needed to establish what sig- 
nature occurred first and therefore where 
the onset begins. 


Why? 

The "where" and "when" questions of sub- 
storm triggers outlined above speak to an un- 
certainty in our empirical view of substorms. 
These are needed before we can assuredly 
seek a clear association with possible exter- 
nal variations. Establishing these morphologi- 
cal characteristics are critical also for 
addressing the more fundamental physical 
questions of what the trigger mechanism is 
and why it operates at a given location and 
time in the substorm cycle. These underlying 
questions should become clearer as the phe- 
nomenology is clarified in the ISTP era. The 
substorm is one of the two most prominent 
dynamic responses of the terrestrial space 
physics environment to solar-induced vari- 
ations (the other being magnetic storms) and 
is the most frequent one. Substorms are an 
important element of "space weather" and 
thus has not only a fundamental research mo- 
tivation, but also a potentially practical appli- 
cation to space- and ground-based systems 
that may be vulnerable to variations of the 
natural space environment. By any measure, 
physical models are only as good as their abil- 
ity to accurately predict a given phenome- 
non. At this stage, we do not understand the 
physics governing the substorm trigger well 
enough. This research area will undoubtedly 
yield the next generation of substorm models 
that go well beyond the ability to describe 
substorm phenomenology, but are truly pre- 
dictive. 
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We combine two theoretical models of plasma sheet particle transport to trace plasma 
sheet particles from the nightside magnetospheric equatorial plane to low altitudes. We 
predict that the low-altitude signature of the plasma sheet is manifestly different at different 
local times. The Guiding Center Transport Model (GCTM) (Onsager el al , 1993) uses 
a given particle energy, position, and an assumed magnetic and electric field to trace the 
low-altitude particle s trajectory to the plasma sheet. The Finite Tail Width Convection 
Model (FTWCM) (Spence and Kivelson, 1993) uses two Maxwellian plasma sources, one 
down! ail and one on the dawnside boundary layer, to generate distribution functions across 
the tail plasma sheet. Using Liouville’s theorem, the distribution functions at the equatorial 
level are projected to low altitudes. When this signature of the plasma sheet is projected 
to low altitudes near the Earth, the spatial variation in the equatorial plasma properties 
produce a latitudinal variation in the isoflux contours. Isoflux contours at high latitudes at 
dawn extend to lower latitudes at dusk, reflecting a larger plasma sheet plasma pressure at 
dusk than at dawn. This prediction is based solely on a. projection of spatial variations in 
equatorial plasma properties. 

1 . Introduction 

The macroscopic structure of the terrestrial plasma sheet is central to many large-scale as- 
pects of magnetotail physics. For instance, in simple magnetohydrodynamic descriptions of the 
quasi-static tail, spatial gradients of the plasma sheet plasma pressure at the magnetic equatorial 
plane are a proxy for tail lobe magnetic field strengths (e.g., Spence et al, 1989). Another example 
is time-dependent magnet ot, ail convection models; the distributions of number density and pres- 
sure constrain magnetic-flux-tube volumes and hence magnetospheric convection (e.g., Erickson, 
1992). In such descriptions of the plasma sheet, plasmas and magnetic fields are intimately linked 
and each provides unique yet coupled insights on the physics of the magnetotail. 

Both theoretical models and empirical models may be used to describe the tail. Empirical 
models of the magnetic field in the magnetotail have improved steadily over the past two decades 
and are now used commonly as a research tool (Peredo et al, 1993). Similar models of magnet otail 
plasmas are not in common use for several reasons: incomplete energy coverage of measured 
spectra; limited orbital coverage; and temporal /spatial ambiguity, since the transit time of a 
large apogee satellite is greater than the dynamic timescale, even during quiescent or quasi-static 
periods. 

This last shortcoming is avoided at intermediate to low altitudes (^3 Re), where polar or- 
biting spacecraft traverse the lower altitude extension of the plasma sheet (the so-called bw horns ) 
relatively quickly. These rapid passages yield “snapshots 1 ' of the equatorial plasma sheet. Stud- 
ies have utilized this approach successfully (Suszcynsky et al . , 1993; Onsager and Mukai, 1995). 
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Despite the advantage gained by low altitude measurements of the plasma sheet,, a single low- 
altitude spacecraft provides only one passage through the plasma sheet horns per orbit; and it 
occurs at only one fixed local time (equivalent to a single equatorial trajectory through the mag- 
netotail plasma sheet). To construct an empirical two-dimensional, instantaneous map of plasmas 
in the equatorial plane requires many satellites traversing the plasma sheet horns simultaneously 
at various local times. Such multi-satellite studies have not been available routinely. 

In this paper, we complement single spacecraft low-altitude data studies. We combine two 
published theoretical models of magnetospheric particle transport to predict the low-altitude sig- 
nature of the plasma sheet as a function of local time. One model specifies the plasma thermal 
pressure and density throughout the equatorial magnetotail. The other specifies the mapping 
of plasma properties between the equatorial plane and intermediate altitudes. The magnetotail 
plasma model is quasi-static yet predicts spatial gradients of the macroscopic plasma quanti- 
ties. These cross- and down-tail variations at the equatorial plane produce differences in the 
low-altitude plasma sheet signatures at different local times. It is important to stress that the 
predicted local time differences arise from spatial differences, a conclusion not easily reached from 
a single spacecraft passage through the plasma sheet horns. We suggest that model syntheses of 
this type improve the interpretation of low-altitude plasma sheet signatures. 

2. Models 

2. 1 Guiding center transport model 

The Guiding Center Transport Model (Onsager et al, 1993) traces guiding centers of ions and 
electrons between low altitudes and the magnetic equatorial plane. The particle guiding centers 
are traced backwards in time to determine the source locations in the equatorial plane of particles 
detectable by a low-altitude spacecraft. For the purposes of this research, we have considered only 
protons. Protons are followed as they travel through specified electric and magnetic fields; their 
energy and magnetic moments are conserved as their guiding centers are traced. The electric 
field is an imposed dawn to dusk field across the tail; magnetic field lines are assumed to be 
equipotentials. The magnetic field model used is the Tsyganenko model (Tsyganenko, 1987). 
The magnetic field model is parameterized by the Kp index and the dipole tilt angle. 

In Fig. 1, we show a typical guiding center trace in the Xqsm-Zgsm plane. The dashed 
line is the trace of the guiding center. The solid lines are the magnetic field lines. The particle s 
guiding center trajectory in two dimensions is determined by the integral effects of its E x B 
drift velocity across magnetic field lines and its velocity along the magnetic field lines. The 
particle being traced “backwards” from an altitude of 10,000 km and at ~73° magnetic latitude 
is a 2 keV proton with a 10° pitch angle. These specify the initial velocity components and 
following the guiding center to its origin, we find that this particle originated in the plasma sheet 
at Xqsm = -53 Rp- In this example of how the GCTM maps guiding centers, the dipole tilt 
angle was —29° and the Kp — 3. 

The inputs to the GCTM are the position, pitch angle (here assumed to be 10°), and energy 
of the particles detectable at low altitudes. The positions are chosen to represent a polar orbiting 
satellite trajectory at several fixed magnetic local times. For the results described here, we 
have assumed an electric field that corresponds to an E x B drift speed of 1 km/s equatorward 
at the initial low altitude position (~ 10.000 km). This electric field when mapped along the 
assumed equipotential magnetic field to the equatorial plane yields an electric field of about 0.08 
mV/ni at Xqsm — —69 Re and about, 0.15 mV/m at X G sm — —15 Re- These values of the 
electric field correspond to earthward convection speeds in the plasma sheet of about 50 km/sec 
at Xqsm — —60 Re and about 25 km/sec at Xqsm = — 15 Re- The inputs to the magnetic 
field model are the dipole tilt angle and the A> used for Fig. 1. 

The model specifies the particle’s position in the equatorial plane and its energy, but not its 
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Fig. 1. A cartoon of the Guiding Center Transport Model. The dashed line represents a 2 keV proton originating 
in the plasma sheet at X = —53 Re detected at ~10,000 km in altitude at 73° latitude. The K v — 3; the 
Earth’s dipole is tilted —29°. 



Fig. 2. A cartoon of the Finite Tail Width Convection Model. There are two sources of plasma, both Maxwellian: 
the dawnside low-latitude boundary layer and the tail source. Particles drift from the source regions and the 
two sets of particles at a given point are then summed (Spence and Kivelson, 1993). 
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flux. To determine the proton flux we must know the plmse space density at the equatorial plane. 
The plasma conditions throughout the plasma sheet equatorial plane' have been calculated using 
the Finite Tail Width Convection Model (FTWCM) (Spence and Kivelson, 1993) (hereinafter 
referred to as SK93). 

2.2 The finite tail width convection model 

The FTWCM is a bounce averaged drift model for ions. It describes the bulk plasma proper- 
ties in the magnetotail from —10 < XasM < —20 and | Ycsm I < 20 R // at the magnetic equator. 
This model has been described extensively elsewhere (Kivelson and Spence, 1988; Spence and 
Kivelson, 1990, and SK93) and is summarized only briefly here. 

The model has two sources of particles which populate' the plasma sheet: the low-latitude 
boundary layer and the distant magnetotail. Figure 2 is a cartoon depicting the model and its 
sources. The model has five input parameters: proton number density and temperature for both 
plasma sources and the cross-tail electrostatic potential. Both plasma source distributions are 
assumed isotropic and Maxwellian. An additional parameter that the model requires is the flux- 
tube volume along the midnight meridian. This is calculated using the Tsyganenko 198^ long 
model (Tsyganenko, 1987). 

Particles drift both earthward and from dawn-to-dusk under the influence of E x B and VB 
drifts from the source distributions. For a finite width tail, only those protons starting in the deep 
tail (dawn boundary layer) with velocities less (greater) than some limiting velocity can reach 
any given final position. At this final location, the ensemble distribution function is the sum of 
those portions of the source distributions that can gain access. The plasma moments have' finite 
velocity limits and this modifies densities and pressures as a function of position within the tail. 

Figure 3 illustrates the model outputs: plasma density and thermal pressure calculated as 
a function of position in the magnetic equatorial plane. For our study, representative quiet 
magnetospheric conditions were used. The model inputs are as follows: a cross tail potential of 
15 kV; a far-tail plasma source temperature, kT T = 2.5 keV and number density. n T = 0.1 cm-*; 
and a low-latitude boundary layer temperature, hT L = 0.7 keV and number density. n L = 0.5 
curl These values are identical to those motivated by and used in SK93. The top panel of 
Fig. 3 shows number density contours in number/cnr* plotted in the equatorial plane, where the 
vertical axis is Yqsm and the horizontal axis is Xqsm ■ The middle panel shows thermal pressure 
in iiDynes/cm 2 . The bottom panel shows the temperature contours in keV. 

Note that each panel demonstrates the strong ^-dependent variation owing to adiabatic 
compression as well as a {/-dependent variation owing to gradient drift effects. The earthward 
increases in density and pressure are well studied and understood in terms of convective transport 
and stress balance (Lui et al , 1994; Spence et a/., 1989). The crosstail variation has been less 
studied empirically and warrants discussion in the context of this paper. 

The large crosstail plasma pressure gradient must be in stress balance with the coincident 
magnetic geometry in the magnetostatic limit. This is not true for the FTWCM using the Fsy- 
ganenko 1987 long model, as the magnetic field was not self-consist. ent.ly readjusted according to 
the predicted plasma pressures. Nevertheless, we believe that the FTWCM provides a reasonable 
first-order approximation. For example, the model cross-tail pressure gradient was considered 
by SK93. They reasoned that pressure balance cannot lx; maintained by an equal and opposite 
crosstail lobe magnetic pressure gradient. Instead, they proposed a magnetic pressure gradient 
within the plasma sheet which increases from dusk to dawn. Evidence for this condition was 
found not only in tail plasma sheet plasma pressure surveys (Liu and Rostoker, 1991), but also in 
magnetic surveys (Fairfield, 1986). The empirical models compare favorably with the FTWCM 
model results (see SK93). We appeal to this explanation for consideration of the appearance of 
this cross- (and down-) tail variation at low altitudes. 
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3. The Melding of the Models 

The GCTM output gives a position and an energy in the near-tail equatorial plane corre- 
sponding to the source location of a particle detected along the low-altitude spacecraft trajectory. 
At these positions, the FTWCM can then be used to determine the distribution function as a 
function of energy. Using Liouville’s theorem, the values of phase space densities in the equatorial 
plane can be assigned to the corresponding locations at low altitudes for comparison with a polar 
orbiting spacecraft, such as Akebono. Such comparisons have been made recently by Onsager 
and Mukai (1995). In their work, they arbitrarily specified phase space densities at the equatorial 
plane and restricted their attention to one magnetic meridian (local time) at low altitudes. In 
this paper, we extend their approach by investigating the low-altitude plasma sheet signature at 
several local times in the context of a two-dimensional physically- motivated model of plasma at 
the equator. 

In combining the models, the assumptions and inputs are as consistent as possible. For 
example, the magnetic field model used is the same for both models. Other aspects of the 
models differ slightly in approach but these differences do not affect the results of the mappings 
qualitatively. For example, the electric fields are treated somewhat differently in each model. 
The electric field magnitudes at the equatorial plane are similar at the down-tail boundary of 
the FTWCM, but differ by a factor of two at the earthward boundary. We stress that for the 
GCTM results presented here, the calculated low-altitude spectra are governed primarily by the 
magnetic field mapping between the low-altitude locations and the plasma sheet. In these cases, 
the details of how the electric field is distributed are a secondary effect and are not treated in 
this study. 

To illustrate the model synthesis we chose to generate synthetic spectra that a low-altitude 
spacecraft plasma analyzer (such as on the Akebono spacecraft (Mukai et al , 1990)) would be 
able to observe. We chose Akebono because it has been used recently in such studies (Onsager 
and Mukai, 1995). Accordingly, the energy bins used for the GCTM were chosen to match 
Akebono energy bins and initial low-altitude locations were chosen along orbital tracks of a polar 
orbiting spacecraft at three local times. The FTWCM parameters were chosen to match the SK93 
values. The synthesized model results are thus representative of relatively quiet magnetospheric 
conditions. The results shown in Fig. 4 demonstrate the prediction of the plasma sheet’s projection 
to various low-altitude locations in a general sense. Of course, any comparison of the model with 
real data would require selection of model parameters consistent with actual conditions. 

Figure 4 shows the combined model output of spectra that we predict would be measured 
by three spacecraft simultaneously traversing the “horns” at three local times: ~2200 hours, 
midnight, and ~0200 hours. The vertical axis is energy measured on a logarithmic scale. The 
horizontal axis is magnetic latitude. The contours are differential energy flux. White on this scale 
represents flux levels below the simulated detector threshold here consistent with the Akebono 
Low Energy Particle (LEP) sensors. 

Figure 4 describes the predicted signature of the plasma sheet at low altitudes as a function 
of local time. This figure illustrates that low-altitude data can be used to remotely sense plasma 
sheet conditions. For example, at any particular local time pass, the flux intensity increases 
as the spacecraft moves to lower latitude reflecting the earthward compression of plasma and 
the commensurate rise in number density. Also, the energy width of the flux contours increases 
equator ward, reflecting the earthward increase of plasma sheet pressure and temperature. The 
flux levels can be used quantitatively to infer plasma sheet bulk properties. 

Of more interest for this paper, though, are the predicted local time dependences which also 
mirror plasma sheet conditions. The important feature on these plots is the apparent latitudinal 
position of the outer edge of the plasma sheet signature (identified by the first observable isoflux 
level) at the three local times. The y-dependence of the plasma sheet temperature (see bottom 
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Fig. 4. The model output compared at three local times: midnight , 0200, and 2200. White represents differential 
energy fluxes below the threshold. 


panel of Fig. 4) yields a clear local time dependence at low altitudes: the flux energy width (a 
measure of temperature) increases from the top to bottom panel of Fig. 4, consistent with the 
dawn-to-dusk energy gradient . This is most clearly seen when comparing the profiles at 2200 and 
0200 magnetic local time where the magnetic mappings, which affect the latitudinal appearance 
of the plasma sheet, are identical within the symmetric Tsvganenko model. The mapping of the 
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Tsyganenko model also explains why the midnight mapping of the first isoflux contour is at lower 
latitudes, since the field is more stretched near midnight. 

4. Conclusions 

Two models were linked: the finite tail width convection model and the guiding center 
transport model. Their output makes specific predictions for plasma sheet ion fluxes detectable 
at low-altitudes. Specifically, before local midnight, the plasma sheet appears to be hotter and 
the particle isofluxes are visible at higher invariant latitudes than after local midnight. 

The spectra in Fig. 4 show the spatial variations in the plasma sheet, signature; the invariant 
latitude where the protons have the same flux differ with local time as does the energy width of the 
measured flux. The importance of this model is that it is time-independent. Therefore, all model 
predictions for differences in the plasma sheet signature are due to spatial effects: the projection 
to low altitudes of the plasma pressure asymmetry in the equatorial plane. This asymmetry is 
due to the differential drift of plasma across the tail. This type of variation should be readily 
apparent in concurrent satellite passes at low or middle altitudes at several loc al time s. 

These models can be used to compare near-simultaneous measurements by multiple space- 
craft. Currently we are in the process of assessing the models against ISEE-2 data (Hirsch et 
a l 1995). The ideal test case is simultaneous measurements at different local times from two or 
more spacecraft. During the International Solar-Terrestrial Physics program it is of great impor- 
tance to have diagnostic tools such as this one to develop links between low and middle altitude 
crossings of the plasma sheet horns and the equatorial plane. 
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